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Summary. — Recent measurements of the correlation between the spin of the 
top and the spin of the anti-top quark produced in proton anti-proton scattering 
at a centre of mass energy of -/s — 1.96 TeV by the CDF and DO collaborations 
are discussed. Using up to 4.3 fb~^ of data taken with the CDF and DO detectors 
the spin correlation parameter C, the degree to which the spins are correlated, is 
measured in dileptonic and semileptonic final states. The measurements are found 
to be in agreement with Standard Model predictions. 



1. Introduction 

Top quark physics at liadron colliders plays an important role in testing the Standard 
Model of particle physics and its possible extensions. 

In the Standard Model the top quark has a very short lifetime, T1/2 ~ 5 x 10"^^ s, 
therefore the definite spin state in which the top anti-top pair is produced is not spoilt 
by hadronisation effects. As a result, the direction of the spin of the top quark is re- 
flected in the angular distributions of its decay products. In contrast to this, the spin 
of light quarks will flip before they decay, making the spin state they are produced in 
unobservable. Furthermore, the theoretical calculations necessary in order to predict 
the angular distributions can be performed for top pairs, resulting in precise theoretical 
predictions which can be tested by experiment. New physics in either the production 
or decay mechanism would modify these angular distributions, making spin correlations 
sensitive to new physics. 

Until recently only one measurement of spin correlations has been performed. Using 
125 pb"^ of data taken during Run I of the Tevatron collider at Fermilab the DO col- 
laboration measured a correlation coefficient in agreement with the Standard Model [1]. 
However, since the sample contained only six events, the sensitivity was too low to rule 
out the hypothesis of no spin correlations. Recently the CDF and DO collaborations 
performed measurements using up to 4.3 fb"^ of data taken with the CDF and DO [2] 
detectors, the results of which are discussed below. 
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2. Observables 

In strong interactions the top and anti-top quark are produced unpolarised at hadron 
colliders, however the tt system is in a definite spin state. At the Tevatron about 85 %, at 
next to leading order, of top quark pairs are produced via quark anti-quark annihilation. 
At threshold these tt systems will be in a ^Si state, whereas the 15% of top quark 
pairs produced via gluon fusion will be in a ^5*0 state. In the first case the top and 
anti-top quark will tend to have their spins parallel, in the second case they tend to be 
anti-parallel. One therefore expects to observe a correlation between the direction of the 
spins. 

The strength of the correlation due to the production mechanism can be expressed as 
the asymmetry A, 

^^ ^^ _ iVn + Nu - iVit - N^i 

Nn + Nu + Ni^ + Nn 

between the number of events with spins parallel, N-^-f and N^^, and the number of events 
with spins anti-parallel, N-^i and iV|^-|-. 

In order to measure the direction of the spin vector a quantisation axis needs to be 
defined. At the Tevatron three sets of quantisation axes, referred to as "spin basis", are 
commonly used. They are shown in Figure 1. 

The simplest is the so called "beamline basis" in which the direction of one of the 
incoming hadrons is used as quantisation axis. This basis is easy to construct and is 
optimal for tt systems produced at threshold. The production asymmetry has been 
calculated at next to leading order (NLO) in QCD a,s A ~ 0.777 [3]. 

The second basis is the "helicity basis" in which the momentum of the (anti)top quark 
in the top-anti-top quark zero momentum frame is used to quantise the (anti)top quark 
spin. At the Tevatron the strength of the correlation is smaller than in the "beamline 
basis", in NLO QCD A = -0.352. The opposite sign arises due to the fact that the spins 
tend to be anti-parallel in this basis. 

Finally the third basis is the "off-diagonal basis". The direction of the quantisation 
axes are defined by the angle lu with respect to the (anti)top quark momentum. The 
angle w is given by tanw ~ -^/l — /3-^ ta.n9, where /3 is the speed of the top quark and 
9 is its scattering angle. This basis interpolates between the "beamline basis" close to 
threshold (low f3) and the "helicity basis" above threshold (large /?). The production 
asymmetry is A = 0.782. While this is slightly larger than in the "beamline basis" it is 
more complex to reconstruct. 

The angular distribution of decay product i in the top quark rest frame is given by: 

1 d(T 1 , 

2 --; - = -{l-a,-cos0,) 

a a cos Oi 2 

where Oi is the angle between the direction of flight of decay product i and the direction 
of the spin vector; oii is the so-called spin analysing power. From Equation 2 it is clear 
that the angular distribution of a decay product with Ui — contains no information 
about the direction of the top quark spin and the angular distribution of a decay product 
with Q!i = ±1 will contain most information. The spin analysing power of the various top 
quark decay products are listed in Table I. The particles with the highest spin analysing 
power are the lepton and the down type quark from the W boson decay. 
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Figure 1. - The three choices of quantisation axis used at the Tevatron. The "beamhne basis" 
(left) is optimal for top pairs produced at threshold, the "helicity basis" (centre) is used for 
above threshold top pairs and the "helicity basis" (right) interpolates between the two. 



In order to observe a correlation between the direction of the spin of the top and 
anti-top quark one must consider the angle of a decay product of the top quark and the 
angle of a decay product of the anti-top quark simultaneously. The double differential 
distribution for a top quark decay product i and anti-top quark decay product j is given 
by: 



(3) 



1 



dV 



1 



„ , „- = - (I — AaiUj cosOiCosOj) 



where a is the total cross section, A is the production asymmetry, and at^j is the spin 
analysing power of the i, j-th decay product. In all analyses presented here the spin 
correlation parameter C — AaiOj is measured. A measurement of the distribution given 
in Equation 3 should be performed as follows: 

1. Reconstruct the top and anti-top quark momenta in the laboratory frame, 

2. Perform a boost from the laboratory frame to the rest frame of the tt system. 
Define the vectors bi and bj along which to quantise the top and anti-top quark 
spins respectively. 

3. Boost the top (anti-top) quark decay product to the top (anti-top) quark rest frame 
and calculate cosOij — bi^j ■ qi,j. 

The difference between the case of no spin correlations, A = 0, and SM spin correlations 
as measured in the "beamline basis", A — 0.777, using leptons as spin analysers is shown 
in Figure 2. 



Table I. - Spin analysing power of the top quark decay products. The up type quark, down type 
quark, neutrino and lepton are the decay products of the W boson. For the antiparticles the sign 
is reversed. 



lepton, down type quark neutrino up type quark 



b quark 



analysing power a 



-1 



-0.31 



-0.31 



-0.41 
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Figure 2. - The distribution cos Oi cos 62 for a sample of top anti-top quark events using generated 
partons. With spin correlations (red dashed) and without (black solid). Here both top quarks 
decayed to leptons which subsequently were used as spin analysers [4]. 



3. Measurements 

While in theory the down type quark is as powerful a spin analyser as the lepton, it 
is more difficult to identify in practise. This leads to two different approaches. In the 
first one, one selects a pure sample of top pairs in which both the top and anti-top quark 
decay to leptons. In the second a sample with higher statistics is selected by requiring 
only one top quark to decay to a lepton. In the following the advantages, challenges and 
results are discussed for both approaches. 

3'1. Dilepton final states. - The advantages of the dilepton final state are that it is 
simple to identify the final state particles of interest and the high purity of the sample. 
The disadvantage is that one suffers from a low branching ratio and needs to deal with 
two neutrinos when reconstructing the kinematics of the event. Both CDF and DO select 
events with two high pt leptons of opposite charge and at least two jets. The detailed event 
selections are described in References [5, 4] for CDF and DO, respectively. In final states 
with same flavour leptons (e"'"e" and fJ-'^ fJ-') the main background arises from Drell-Yan, 
Z/7* — >■ t£^, production. In the e/u final state the main background is instrumental, 
this occurs mainly due to W+jets events in which a jet is misidentified as a lepton. The 
second largest background arises due to semileptonic decays of Z/7* — >■ t't^. Further 
sources of background in all three final states are the diboson processes WW, WZ and 
ZZ. Both signal and background are modelled using Monte Carlo simulation, except for 
the instrumental background which is estimated from data. 

In order to reconstruct the momentum of the top and anti-top quark, one needs to 
deal with the two neutrinos in the final state. To fully characterise the kinematics of 
the final state one needs 18 quantities, assuming the masses of the final state particles 
are known. While the leptons and jets are observable in the detector, the two neutrinos 
escape detection. It is possible to infer the sum of the momenta of the neutrinos in the x 
and y plane from the missing transverse energy, ^^ and ^|, . Using this information and 
making an assumption about the mass of the W boson and the top quark it is possible 
to write down a set of quartic equations which fully describe the final state. Solving 
them yields up to four solutions per event. Additionally one needs to try both lepton-jet 
pairings which increases the number of possible solutions to eight. 
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In the CDF measurement a likelihood function is constructed from several observables 
and maximised with respect to the unknown neutrino momenta (p^, pp) and the energies 
of the bottom quark jets {Ef'^^'', Ef'^'^y. 

L (p., p„ Er^\ i?f--) ^ P (p*/) P (p**) P (M,f) X 



b 



, , rimcas i-igucss \ 2\ / / pmeas pgucss ^ 2^ 



^_p™ _i(£r_£rr: x^— exn -i 

where P (p'j^) , P (pt) ^^^ P {Mti) are probability density functions obtained from Pythia 
tt Monte Carlo events, E^^^ the measured energies of the bottom/anti-bottom quark 
jets, fi^T^ the measured components of ^t, and ai the respective resolutions. The max- 
imisation is performed for both lepton-jet pairings and the combination with the larger 
L is kept. 

As the Pythia Monte Carlo simulation does not contain spin correlations, templates 
for values of C — —1, 0.8, . . . , 1 are obtained by reweighting the signal Monte Carlo at 
the generator level using a weight w ~ 1 — C • cos 9i cos 62 ■ For each value of C a two 
dimensional template in the decay angles of the lepton and anti-lepton, cos 61+ , cos Og- 
and a template in the decay angles of the bottom quark jets, cos 9b , cos O^ is created. 
The two templates are fit with an analytical function /^'^ (x, y; C). The measurement 
is then performed on the N candidate events by maximising the likelihood function: 

N 



LiC)=l[fix,y:C)f{x,y;C) 



i=0 

In order to extract limits from the measurement, a confidence belt according to the 
Feldman-Cousins prescription [6] is created. This naturally includes both statistical and 
systematic uncertainties and allows one to decide before looking at the data whether 
to quote a one or two sided limit. Using 2.8 fb"^ of data the best fit value is C = 
0.32tQyg(stat + syst) and the corresponding confidence belts are shown in Figure 3. 
The measurement was performed in the "helicity basis". The result is consistent with the 
expected value of C = 0.782. The largest contributions to the systematic uncertainty 
come from evaluating the PDF uncertainties and the finite number of Monte Carlo events 
used to form the templates. 

At DO the neutrino weighting technique is used to solve for the event kinematics. 
By making an assumption about the rapidity, 77, of the neutrino and anti-neutrino, it is 
possible to solve the event kinematics, while not using _^y and fiJj^ in the process but 
instead to assign a weight, w, to each solution given by: 

. . exp (jrr-'^.-^.f] , ,,p (jK-^y-^yf) 



where i^x.y and Dx,y are the x and y components of the neutrino and anti-neutrino 
momentum for a given solution and a is the f]^ resolution. Many solutions are obtained 
by sampling the neutrino and anti-neutrino rapidity based on Monte Carlo simulation. 
No dependence of the neutrino rapidity on the presence of spin correlations is observed. 
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Figure 3. - The 68% (stat only), 68% and 95% Confidence Level intervals constructed according 
to the Feldman-Cousins prescription including statistical and all systematic uncertainties for 
the CDF measurement. The best fit value is C = 0.32l°;^^ [5]. 



The weighted mean of all solutions for an event is used as estimator for the true value of 
cos0£+ cos 9 1- . 

As for the CDF measurement, the Pythia Monte Carlo simulation is used to model 
the signal sample. A one dimensional template in the variable cos 6*^+ cos 6i- is created 
for C = and C = 0.777 by reweighting the distribution at the generator level. In order 
to extract a value of C a linear combination of the two templates is fit to the data. 

Pseudo-experiments are created for each value of C and fit with signal and background 
templates. Each source of systematic uncertainty is considered as a nuisance parameter 
during the fit. Feldman-Cousins confidence belts are constructed from the pseudo exper- 
iments. Using up to 4.2 fb"^ of data the best fit value is C = — 0.17^Q'53(stat + syst). 
In this measurement the "beamline basis" was used and the measured value is consistent 
with the Standard Model expectation of C = 0.777 at the two sigma confidence level. 

The two main sources of systematic uncertainty are the variation of the assumed 
top mass during the event reconstruction from 175 GeV to 170 GeV and the test of the 
reweighting method. For the latter, the two Pythia signal templates were replaced by 
Alpgen, which contains spin correlations, and MC@NLO where spin correlations were 
turned off. 

3'2. Semileptonic final states. - Selecting semileptonic events results in a higher yield, 
but the challenge is to identify the down type quark. This is done probabilistically by 
choosing the jet closest to the bottom type jet in the W boson rest frame [7], which will 
result in picking the correct jet about 60% of the time. 

Events are selected by requiring at least one high pT, central lepton, large missing 
transverse energy and four or more jets, one of which must be identified as a b-jet. The 
backgrounds are estimated both from simulation and data. For details of the selection 
see Reference [8]. Using 4.3 fb"^ of data a total of 1001 events are selected of which 786 
are expected to be top pair events. 

When produced in pairs the top and anti-top quark either have the same helicity or 
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Figure 4. - Left: the 68%, 95% and 99% Feldman- Cousins confidence belts are sliown. The 
best fit value can be read of at the intersection of the dashed black line and the thin blue line. 
Right: The sum of all dilepton channels is shown. The open black histogram shows the expected 
distribution for the case of no spin correlations, C — and the filled red histogram the expected 
distribution for Standard Model spin correlations, C = 0.777 [4]. 



opposite helicity. The fraction of top pairs with opposite helicity is given by: 



fo = 



cr {tlitji + ti^tL + t^tL + tLtji) 



where u {tL.RtL,R) denotes the cross section for each possible helicity configuration. 
Using Equation 1 one can show that a measurement of fo is equivalent to a measurement 
of A in the helicity basis. 

One template for top pairs with same helicity and one template for top pairs of 
opposite helicity are created using a modified version of the Herwig event generator. 
The opposite helicity fraction is extracted with a binned maximum likelihood fit of the 
two templates to the data, with contributions from backgrounds taken into account. 
The best fit value is fo = 0.80 ± 0.26(stat + syst) or equivalently A = 2fo — 1 = 
0.60 ± 0.52(stat + syst). This is consistent with the Standard Model expectation of 
A = 0.4. The two main systematic uncertainties are Monte Carlo statistics and jet 
energy scale. 

4. Conclusions 



The spin correlation parameter C has been measured in dilepton and semileptonic 
decays of top and anti-top quark pairs using up to 4.3 fb"'^ of data collected with the 
CDF and DO detectors. Measurements were performed in the "beamline", "helicity" and 
"off-diagonal" bases. The measurements are found to be in agreement with the Standard 
Model predictions. All three measurements are still statistically limited. Considering 
that the Tevatron collider has delivered nearly twice as much integrated luminosity since 
the analyses have been performed, updates of all measurements can be expected soon. 
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Figure 5. - The best fit of same iielicity, opposite iielicity and background templates for tlie 
CDF semileptonic decay ciiannel. On tiie (left) the distribution of the product of the decay 
angle of the lepton and the bottom quark. On the (right) the distribution of the product of the 
decay angle of the lepton and the down type quark. The best fit value from a simultaneous fit 
to both distributions is fo = 0.8 ± 0.26(stat + syst) or C = 0.6 ± 0.52(stat + syst) [8]. 
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